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ABSTRACT 

We present optical spectroscopic and photometric observations of supernova (SN) 2008D, associated 
with the luminous X-ray transient 080109, at > 300 days after the explosion (nebular phases). We 
also give flux measurements of emission lines from the H ii region at the site of the SN, and estimates 
of the local metallicity. The brightness of the SN at nebular phases is consistent with the prediction 
of the explosion models with an ejected ^^Ni mass of 0.07 Mq, which explains the light curve at early 
phases. The [O l] line in the nebular spectrum shows a double-peaked profile while the [Ca ii] line 
does not. The double-peaked [O i] profile strongly indicates that SN 2008D is an aspherical explosion. 
The profile can be explained by a torus-like distribution of oxygen viewed from near the plane of the 
torus. We suggest that SN 2008D is a side-viewed, bipolar explosion with a viewing angle of > 50° 
from the polar direction. 

Subject headings: supernovae: general — supernovae: individual (SN 2008D) — nuclear reactions, 
nucleosynthesis, abundances — line: profiles 



1. INTRODUCTION 

On 2008 January 9, a luminous X-ray transient was 
serendipitously discovered in NGC 2770 during the 
follow-up observation of SN 2007uy in the same galaxy 
with the Swift satellite (Berger & Soderberg 2008). An 
optical counterpart was also discovered at the position of 
the transient (Deng & Zhu 2008; Valenti et al. 2008b), 
and it was named supernova (SN) 2008D (Li & Filip- 
penko 2008). 

The total energy emitted at X-ray wavelengths is ^ 
2 X lO'^^ ergs, smaller than that of long gamma-ray bursts 
(GRBs) by a factor of > 1000 (Soderberg et al. 2008). 
The origin of the X-ray emission is being debated. Soder- 
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berg et al. (2008), Chevalier & Fransson (2008), and 
Katz et al. (2009) interpreted the X-ray transient as a 
SN shock breakout. On the other hand, Xu et al. (2008), 
Li (2008) and Mazzah et al. (2008) suggested that the 
transient is the least energetic end of GRBs or X-ray 
flashes. 

SN 2008D is classified as Type lb because of the pres- 
ence of He lines (Soderberg et al. 2008; Mazzali et al. 
2008; Malesani et al. 2009; Modjaz et al. 2008c) while 
SNe associated with GRBs are all Type Ic (without He 
lines). The progenitor star of SN 2008D has a He layer 
prior to the explosion, and the He core mass is estimated 
to be 6-8 Mq (Tanaka et al. 2009). 

In this paper, we present optical spectroscopic and 
photometric observations of SN 2008D at > 300 days af- 
ter the explosion (nebular phases) with the Subaru tele- 
scope equipped with FOCAS (Kashikawa et al. 2002) and 
the Very Large Telescope (VLT) equipped with FORSl 
(Appenzeller et al. 1998). We find that the spectrum of 
SN 2008D shows a double-peaked [O i] emission profile, 
suggesting that SN 2008D is a bipolar explosion viewed 
from near the equatorial direction. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spectroscopy 

A spectroscopic observation of SN 2008D was per- 
formed on 2009 January 6 (texp = 363 days) with the 
Subaru telescope. Hereafter texp denotes the epoch 
in observer's frame measured from the explosion date, 
MJD=54474.56 (Malesani et al. 2009; Modjaz et al. 
2008c). 

Blue (4800-8000 A) and red (6000-9000 A) spectra were 
taken separately. We used a slit of O'.'S width and two 300 
lines mm~^ grisms. The slit was placed with the posi- 
tion angle of 0° (north-south) . Typical seeing during the 
observation was FWHM ~ 1'.'2 measured with star pro- 
files in R band. No filter was used for the blue spectrum, 
while the 058 filter was used for the red spectrum to 
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TABLE 1 

Observed flux of narrow 

emission lines 



Fig. 1. — 40'.'0 x 40"0 section of the i?-band image centered on 
SN 2008D taken witli Subaru telescope at tcxp = 362.9 days. The 
position of the SN is shown with the arrow. Two rectangles show 
the regions where the galaxy spectrum is extracted. 



eliminate second order light. The exposure time is 9000 
and 3300 sec for the blue and red spectra, respectively. 
The wavelength resolution is A/AA ~ 650. 

The observed data were bias-subtracted and flat- 
corrected, and then two-dimensional spectra were ex- 
tracted. Wavelength calibration was performed with a 
Th-Ar lamp. 

Since the SN is located inside a brightness trough 
which surrounds the SN position with a diameter of ~ 4 
arcsec (see Figure 1), simultaneous subtraction of the 
background galaxy and sky lines using the neighboring 
region leads to negative SN flux. Thus, we first extracted 
a one-dimensional galaxy spectrum using the region at 3- 
12 arcsec north and south from the SN position (marked 
with boxes in Figure 1). Then, the galaxy spectrum was 
subtracted using the spatial profile at 6920-6990 A, in 
which the SN docs not show any emission line. 

To maximize S/N of the SN spectra especially at emis- 
sion lines, we extracted the galaxy spectrum using a wide 
region. Since the global color of the galaxy spectrum de- 
pends on the position somewhat, this leads to a slight 
over-subtraction in the extracted SN spectrum in the 
blue (see Figure 4). However, the profile of the SN emis- 
sion lines are not affected by the choice of integrated 
region for the galaxy spectrum. 

The background-subtracted spectra were summed into 
one-dimensional spectra, and then, the fiux was cali- 
brated using the spectrophotometric standard star Fcige 
34 (Oke 1990). The blue and red spectra were combined, 
weighted by the exposure time. Finally, the spectrum is 
scaled with the i?-band magnitude (Section 2.3). 

2.2. Narrow emission lines from the H ii region 

There are very strong narrow emission lines from the 
H II region at the position of the SN. Because of the 
faintness of the SN at nebular phases, the flux of the 
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H/3 A4861 
[O III] A 5007 
[N ii] A 6548 
Ha A 6563 
N ii] a 6584 
S II] A 6717 
S II A 6731 
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1.97 ± 0.61 
2.32 ± 0.66 
33.1 ± 0.9 
7.00± 1.27 
4.69 ± 0.52 
3.18 ± 0.55 
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narrow emission lines can be easily measured. Table 1 
shows measured flux of the emission lines. Extinction is 
not corrected for. 

From the line ratios, wo can estimate gas metallicity 
at the site of SN 2008D. Thone et al. (2009) presented 
metallicity measurements using the data taken soon after 
the explosion. We use N2 and 03N2 indexes to estimate 
the mctaUicity (Pettini & Pagel 2004): N2 = log([N ii] 
A6584/ Ha) and 03N2 = log( ([O ill] A5007/ H/3)/ ([N ll] 
A6584/ Ha) ). Note that these indexes are almost inde- 
pendent on the extinction because of the vicinity of the 
fines ([N ii] and Ha, [O ill] and H/3). Pettini & Pagel 
(2004) give the following relations calibrated with the 
metallicity estimated by Tg method: 12 -I- log(0/H) = 
9.37 + 2.03XN2 + 1.26xN22 + 0.32xN23 and 8.73 - 
0.32xO3N2, for N2 and 03N2 indexes, respectively. 

Using these relations, metallicity is estimated to be 
12+log(0/H) = 8.5 ± 0.2 (N2) and 8.7 ± 0.2 (03N2). 
The error above includes the la dispersion in the relation 
(0.18 and 0.14 for N2 and 03N2, respectively). These are 
consistent with the solar metallicity (12-|-log(0/H)=8.66, 
Asphmd et al. 2004) within the error, but larger than 
those by Thone et al. (2009) by -0.2dex (with the 
same indexes) . The host galaxies of the GRB-associated 
SNe (GRB 030329/SN 2003dh, GRB 031203/SN 20031w, 
XRF 020903, XRF 060218/SN 2006aj) have 12 + 
log(0/H) < 8.1 in 03N2-Te scale (Modjaz et al. 2008a). 
Thus, the metallicity at the site of SN 2008D is higher 
than that in those environments, as pointed by Thone 
et al. (2009). 

2.3. Photometry 

Photometric observations in BVRI were performed 
on 2008 November 22 (t^xp = 317.9 days), 2009 Jan- 
uary 6 (icxp = 362.9 days. Figure 1), and February 19 
{texp = 406.6 days). We detected a point source at the 
SN position in the R band at three epochs and the I 
band at one epoch. A log of the observations is shown 
in Table 2. We performed PSF photometry for the R- 
and /-band images, where a source is detected. For the 
images in which the SN was not detected, limiting mag- 
nitudes were derived using the local sky noise around the 
SN position, and further checked using artificial sources. 

Since the locally strong emission of the H II region 
is present at the SN position, our i?-band photometry 
overestimates the SN brightness. We have corrected for 
this effect by estimating the contribution of the narrow 
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TABLE 2 
Log of Photometric Observations 



Date MJD Epoch ^ ms my mn ^ mj Mboi Telescope 

2008 November 22 54792.5 317.9 - >22.4'^ 23.0to4 (22.7 ± 0.3) 22.4 ± 0.5 -10.0 ±0.8 Subaru 

2009 January 6 54837.5 362.9 >22.8'= >22.2': 23.8t° g (23.2 ± 0.4) >22.1'^ -9.8 ± 1.2 Subaru 
2009 February 19 54881.2 406.6 - >23.6'= ' (23.9 ± 0.4) >22.9'= - VLT 



''Days after the explosion (MJD=54474.56) in the observer's frame. 

''Value in parenthesis is total brightness without correction of the H II region. 

'^5(7 upper limit. 
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Fig. 2. — Reddening corrected, absolute BVRI light curves of 
SN 2008D. The red horizontal line shows estimated brightness of 
the H II region in R band. 



lines using the observed spectrum at t^xp = 363 days 
(mij^gai ~ 24.1 mag, the red horizontal line in Figure 



2). The error in the SN magnitudes at toxp — 317.9 
and 362.9 days does not include the uncertainty caused 
by the contamination from galaxy light. For the data 
at tcxp = 406.6 days, we refrain from measuring the SN 
component because of the large uncertainty in the total 
flux as well as in the galaxy contribution. 

3. LIGHT CURVES UNTIL NEBULAR PHASES 

Figure 2 shows the absolute magnitudes of SN 2008D. 
For the reddening, E{B — V) — 0.65 mag is assumed 
(Mazzah et al. 2008; Modjaz et al. 2008c), and the red- 
dening law of Cardelli et al. (1989) is adopted. For the 
distance to the SN, /n = 32.46 is assumed (Modjaz et al. 
2008c). The data at early phases {tcxp < 100 days) are 
taken from T. Minezaki et al. (in preparation). 

We estimate the bolometric luminosity of SN 2008D 
at nebular phases. Since multicolor observations are not 
available at nebular phases, we use the bolometric correc- 
tion of the well-observed SN 2002ap (Yoshii et al. 2003; 
Tomita et al. 2006). We define the bolometric correction 
as BC = A/boi — Mr,cov, where Mboi is the bolometric 
magnitude constructed from UBVRIJHK magnitudes, 
and Mji^coT is the reddening corrected, absolute i?-band 
magnitude (as plotted in Figure 2). The BC of SN 2008D 
is found to be consistent with that of SN 2002ap within 
0.3 mag at early phases {BC ~ 0.2 - 0.7 mag). The BC 
of SN 2002ap at nebular phases are apphed to SN 2008D 
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Fig. 3.— Pseudo-bolometric Ught curve (black) of SN 2008D. The 
bolometric magnitudes are derived from U BVRIJHK magnitudes 
at early phases (toxp < 100 days, T. Minezaki, in preparation) 
while they are estimated using the bolometric correction for SN 
2002ap at late phases (Tomita et al. 2006). The solid lines in 
blue and green show the results of radiative transfer calculation by 
Tanaka et al. (2009). The gray line shows the total energy release 
from 0.07 Mq of ^SNi. 



by interpolating the epoch. The estimated BC for SN 
2008D is 0.95 ± 0.08 and 0.75 ± 0.2 mag at t^xp = 317.9 
and 362.9 days, respectively. For the bolometric mag- 
nitudes, a systematic uncertainty of 0.3 mag has been 
added. 

The bolometric light curve is shown in Figure 3. The 
models by Tanaka et al. (2009, blue and green lines), with 
an ejected ^^Ni mass of 0.07 Mq, are shown for compar- 
ison. These models reproduce the bolometric light curve 
at early phases. The bolometric magnitudes at nebular 
phases are consistent with the prediction of these models. 
Ejecta mass and kinetic energy are different in two mod- 
els [(Mej/Mg, ^JK/lO^^erg) = (4.4, 3.7) and (6.2, 8.4) 
for model HE6 and HE8, respectively]. However, since 
the difference in the predicted luminosity is very small, 
our observations at nebular phases do not discriminate 
between these two models. 

4. NEBULAR SPECTRUM 

Figure 4 shows the nebular spectrum of SN 2008D 
(blue line) taken at ^oxp — 363 days after the explosion. 
The spectrum is binned into 10 A, similar to the wave- 
length resolution at 6300 A. It is compared with nebular 
spectrum of SNe 1998bw (broad Ic, Patat et al. 2001), 
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Fig. 4. — Nebular spectrum of SN 2008D (blue line) compared 
with that of SNe 1998bw (broad Ic, Patat et al. 2001), 2002ap 
(broad Ic, Mazzali et al. 2007), 2005kl (Ic, Maeda et al. 2008) and 
2004gq (lb, Maeda et al. 2008). The epochs are given as days since 
the explosion (days since the discovery are given for the two epochs 
marked with * as light curves are not available). The spectra are 
normalized to the peak flux of the [O l] AA6300, 6364 line of SN 
2008D. The spectra are shifted by 16, 12, 8, 4, and from top to 
bottom for clarity. 



2002ap (broad Ic, Mazzali et al. 2007), 2005kl (Ic, Maeda 
et al. 2008), and 2004gq (lb, Maeda et al. 2008). The 
spectrum of SN 2008D clearly shows emission lines of 
[O i]AA 6300,6364 and [Ca ii]AA 7291,7323. These lines 
are commonly seen in the nebular spectra of other Type 
Ib/c SNe. The Ca ii IR triplet is marginally detected 
around 8600 A. Narrow lines at 4800-5000 A and 6500- 
6900 A originate in the H ii region of the host galaxy. 

4.1. Line profile 

The [O i] Hne of SN 2008D shows a double-peaked pro- 
file. However, given the strong narrow emission lines 
from the H II region, the broad [O i] line of the SN could 
possibly be contaminated by the narrow [O i] line from 
the H II region. Since the strength of the narrow [O l] line 
is ~ — 15 of Ha around the SN position, we subtract 
the scaled, gaussian-fittcd Ha profile at the position of 
[O i]. Figure 5 shows the original profile (no correction, 
top) compared with the profiles corrected with 1/30 and 
1/15 of the Ha line {middle and bottom, respectively). 
Although the peak at u ~ km s~^ in the original spec- 
trum may be due to the contamination of the narrow 
emission line, the broad, double-peak profile is not af- 
fected by the contamination. 

The [O l] line shows a double-peaked profile, while such 
clear two peaks are not seen in the [Ca ii] line. Although 
these properties were also seen in the spectrum taken at 
iexp = 109 days (Modjaz et al. 2008c), transparency of 
the ejecta was not sure at such a transition epoch from 
photospheric to nebular phases. In fact, evolution of the 
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Fig. 5.— The [O l] line profile of SN 2008D and possible con- 
tamination by the narrow [O l] line from the H II region. Top: line 
profile extracted as in Section 2.1 without correction of a possible 
narrow [O l] line from the H II region. Middle and bottom: line 
profiles corrected for a possible narrow [O l] line. The strength of 
the narrow [O l] line is assumed to be 1/30 and 1/15 of the Ho 
line, respectively. The gray lines show unbinned spectra while the 
black lines show the spectra binned into 10 A. 



profile from toxp = 109 days is not significant. Possible 
change can be seen in the redder peak of the [O i] line. At 
texp = 109 days it is clearly located a,t v ^ km s""'^, but 
the position can be redder at ioxp — 363 days, depending 
on the contamination of narrow [O l] line from the H II 
region. 

The [O i] line profile is compared with other SNe in 
the left panel of Figure 6. Double-peaked profiles are 
also seen in other Type Ib/c SNe (SoUerman et al. 1998; 
Mazzah et al. 2005; Valenti et al. 2008a; Maeda et al. 
2008; Modjaz et al. 2008b). The fraction of SNe showing 
a double-peaked [O i] line is roughly 40 ± 10% (Maeda 
et al. 2008). 

Since our observation is late enough for the ejecta to be 
optically thin, the double-peaked profile of the [O i] line 
is unlikely to be caused by optical depth effects (Tauben- 
berger et al. 2009). It is also unlikely to be caused by 
a combination of the two [O i] lines (AA 6300 and 6364) 
since the strength of the two peaks is comparable (the 
strength ratio is 3:1 in optically thin limit, Leibundgut 
et al. 1991). 

In addition, if the double-peaked [O i] profile were 
caused by asphericity in the distribution of the heating 
source C^^Ni), the [Ca ii] line would also show a double- 
peaked profile. Thus, the profile of [O i] line reflects the 
distribution of excited O i. The [O i] and [Ca ii] lines may 
arise from different sites (Fransson & Chevalier 1989). 

4.2. Comparison with the model 

A simple explanation for the double-peaked [O i] pro- 
file is that oxygen has a torus-like distribution and our 
line of sight is near the plane of the torus (Maeda et al. 
2002; Mazzali et al. 2005; Maeda et al. 2008; Modjaz 
et al. 2008b, but see also discussion by Milisavljevic et al. 
2009). The torus- like distribution of the presupernova el- 
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Fig. 6.— Left: Profiles of the [O i] AA6300, 6364 lines plotted 
versus Doppler velocity measured from 6300 A. The spectra are 
the same as Figure 4, but SN 2003jd (broad Ic, 330 days Mazzali 
et al. 2005; Valenti et al. 2008a) is added. SNe 1998bw, 2004gq, 
and 2002ap show single-peaked profile while SNe 2008D and 2005kl, 
and 2003jd show double-peaked profile. The gray line for SN 2008D 
shows the spectrum corrected for the possible narrow [O l] contam- 
ination (Figure 5). The spectra are shifted by 17.5, 14.0, 10.5, 7.0, 
3.5, and 0.0 from top to bottom. Narrow lines from the host galax- 
ies are removed. Right: Profiles calculated with model BP8 by 
Maeda et al. (2006, 2008). The number represents the viewing an- 
gle measured from the polar direction. The spectra are shifted by 
16, 12, 8, 4, and from top to bottom. 



ements, such as oxygen and carbon, can be realized in a 
bipolar explosion (see Figure 7 for a schematic picture), 
where more energy is deposited in the polar region, and 
nucleosynthesis takes place there (e.g., Nagataki et al. 
1997; Khokhlov et al. 1999; Maeda & Nomoto 2003). A 
bipolar explosion has also been suggested for SN 2008D 
by polarimetric observations (Gorosabel et al. 2008). 

Nebular emission profiles from bipolar models have 
been examined by Maeda et al. (2006), and compared 
to nebular spectra of various SNe Ib/c (Maeda et al. 
2008). In the right-hand panel of Figure 6, the profiles 
calculated with a bipolar explosion model BP8 (highly 
aspherical model) by Maeda et al. (2006) are shown for 
comparison. The kinetic energy of the model is 2 x 10^^ 
erg (/.„ = 0.7 defined in Maeda et al. 2006), which is 
smaller than that estimated for SN 2008D (e.g., Tanaka 
et al. 2009). This is because the models of Maeda et al. 
(2006) include only the C-^O core. In SN 2008D, and 
in the models of Tanaka et al. (2009), the He layer is 
present, and it contains the rest of the kinetic energy. 

The line profile of SN 2008D is explained most nicely 
by a viewing angle of ^ 60 — 70°. The double-peaked 
profile is seen only when the viewing angle is > 50°. If 
the asphericity is smaller, this dividing angle is larger. 
Since the viewing angle and degree of asphericity have a 
degenerate effect, we cannot argue that a large aspheric- 
ity as in model BP8 is preferred for SN 2008D. However, 
it is true irrespective of the models that our line of sight 
should be > 50° from the polar direction. If the angle 
is smaller, the [O i] line would show a single-peaked pro- 
file (see SNe 1998bw, 2004gq and 2002ap and models in 
Figure 6). 



'heavy elements 
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Fig. 7. — Schematic illustration of a bipolar explosion model 
for SN 2008D. O-rich material forms a torus-like distribution. The 
emssion line from the torus-like distribution has a double-peaked 
profile if it is seen from near the equatorial direction (> 50° from 
the pole). 



5. CONCLUSIONS 

We have presented spectroscopic and photometric ob- 
servations of SN 2008D at nebular phases. Flux measure- 
ments of the narrow emission lines from the H ii region 
at the site of the SN and the estimates of the local metal- 
licity were also given. In photometric observations, the 
SN was detected in R and /. The brightness at nebular 
phases is consistent with the prediction of explosion mod- 
els with an ejected ^^Ni mass of 0.07 Mq, which explain 
the light curve at early phases. 

In the nebular spectrum taken at a sufficiently late 
phase, the [O i] line clearly shows a double-peaked pro- 
file, while such a profile is not seen in the [Ca ii] line. 
The evolution of the profile from tcxp = 109 days (Mod- 
jaz et al. 2008c) turned out not to be significant. The 
double-peaked profile cannot be explained by a spheri- 
cally symmetric explosion, and it strongly indicates that 
SN 2008D is an aspherical explosion. The profile can 
be explained by a torus-like distribution of excited O i 
viewed from the side (Figure 7). Our line of sight is 
> 50° from the polar direction, irrespective of the degree 
of asphericity. 
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